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ABSTRACT 

We study the star formation rate (SFR) - stellar mass (M*) relation in a self-consistent manner 
from < z < 2.5 with a sample of galaxies selected from the NEWFIRM Medium-Band Survey. 
We find a significant non- linear slope of the relation, SFR oc M° 6 , and a constant observed scatter 
of 0.34 dcx, independent of redshift and M*. However, if we s elect only b l ue ga laxies we find a 
linear relation SFR cx M*, similar to previous results at z — by iPeng et al.l (|2010f ). This selection 
excludes red, dusty, star-forming galaxies with higher masses, which brings down the slope. By 
selecting on Ltr/Ltjv (a proxy for dust obscuration) and the rest-frame U—V colors, we show that 
star-forming galaxies fall in three distinct regions of the log(SFR)-log(M*) plane: 1) actively star- 
forming galaxies with "normal" dust obscuration and associated colors (54% for log(M*) > 10 at 
1 < z < 1.5), 2) red star-forming galaxies with low levels of dust obscuration and low specific SFRs 
(11%), and 3) dusty, blue star-forming galaxies with high specific SFRs (7%). The remaining 28% 
comprises quiescent galaxies. Galaxies on the "normal" star formation sequence show strong trends of 
increasing dust attenuation with stellar mass and a decreasing specific SFR, with an observed scatter 
of 0.25 dex (0.17 dex intrinsic scatter). The dusty, blue galaxies reside in the upper envelope of the 
star formation sequence with remarkably similar spectral shapes at all masses, suggesting that the 
same physical process is dominating the stellar light. The red, low-dust star-forming galaxies may be 
in the process of shutting off and migrating to the quiescent population. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Galaxies show a strong correlation between their 
star formation rate (SFR) and stellar mass (M*) 
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galaxies on this "star formation sequence" were forming 
stars at much higher rates in the d istant universe relative 
to today (e.g. JMadau et al1ll996D ; for a given mass, the 
SFR has been decreasing a t a steady rate by a factor of 
- 30 from z ~ 2 to z = dDaddi et alj 120071) . although 
it appears to be roug hly constant from z ~ 7 to z ~ 2 
([Gonzalez et al.l [20101 ). The star formation sequence is 
obser yed to have a rough ly constant scatter out to z <~ 1 
(e.g., iNoeske et afll2007l) . 

Generally, star formation is thought to be regulated 
by the balance between the rate at which cold gas is ac- 
creted onto the galaxy an d feedback (e.g.. iDutton et all 
I2010L iBouche et ahll2010( ). whereas the relation between 
gas surface density and SFR surface density does not 
appear to change, at least out to th e highest redshifts 
accessible to molec ular gas studies (jDaddi et al.1 120101 : 
iTacconi et al.l l2010l) . The star formation sequence may 
be a natural consequenc e of 'cold mode accretion' (e.g., 
iBirnboim fc Dekel 120031 ). as the SFR is approximately 
a steady function of time and yields a relatively tight 
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relationship between SFR and M*. Feedback may ef- 
fect the slope of the star formation sequence, whereas 
the evolution of the normalization is thought to result 
from evolution in gas densities with redshift. The scat- 
ter in the star formation sequence may reflect variations 
in the gas accretion history, and is predicted to be insen- 
sitive to stellar ma ss, redshift and feedback efficiencies 
(jDutton et al.|[2010h . 

With a wealth of multi-wavelength data, includ- 
ing ultraviolet (UV) to near-infrared (NIR) photo- 
metric surveys, Ha spectroscopic surveys, and mid- 
to far-infrared photometry from Spitzer /MIPS at 
24/im,.ffersc/ieZ/PACS 70-160/Ltm, and Herschel /SPIRE 
250-500^m, much work has been done to calibrate 
differ ent SFR indicators over broad redshift ranges 
(e.g-JMuzzin et al.l20"l0tlHwang et al.l2010UWuvts et all 
I2011allbl : IReddv et al.l I2012D . Uncertainties in the star 
formation sequence may now be dominated by selec- 
tion effects and observational biases. Selection effects 
can be important, as it is well known that a signifi- 
cant fraction of galaxies have very low SFRs, well below 
the star formation sequence. Furthermore, there exists 
a population of dusty star-forming galaxies with sim- 
ilar rest- frame optical colors t o these quiescent galax- 
ies fe.g.. iBrammer et al.l l2009h . The selection of star- 
forming galaxies and the treatment of quiescent galaxies 
can influence the measured relation and its scatter. 

With the accurate photometric redshifts and photom- 
etry of the NEWFIR M Medium-Band Survey (NMBS; 
I Whitaker etlIll20Tl . we can now extend studies of the 
star formation sequence out to z = 2.5 for the largest 
mass-complete sample of galaxies to date. In this letter, 
we study the evolution of the log(SFR)-log(M*) relation 
from < z < 2.5, selecting star- forming galaxies in a 
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Fig. 1. — The SFR-mass sequence for star-forming galaxies has a non-linear slope out to z = 2.5 (dotted line is linear). The running 
medians and scatter are color-coded by redshift, with a power law fit above the mass and SFR completeness limits (solid lines in bottom, 
right panel). 



way that is independent from the SFR indicator. We 
show that star-forming galaxies with different colors and 
Lir /Ltjv ratios occupy different regions of the log(SFR)- 
log(M*) plane. 

We assume a ACDM cosmology with fiM=0.3, 
S1a=0.7, and Hq=70 km s" 1 Mpc^ 1 throughout the pa- 
per. 

2. DATA AND SAMPLE SELECTION 

Our sample of gal axies is drawn from the NMBS 
(jWhitaker et alJ[20Ll . This survey employs a new tech- 
nique of using five medium-bandwidth NIR filters to sam- 
ple the Balmer/4000A break from 1.5 < z < 3.5 at 
a higher resolution than the standard broadband NIR 
filters. The combination of the medium-band NIR im- 
ages with deep optical medium and broadband photom- 
etry and IRAC imaging over 0.4 deg 2 results in accu- 
rate photometric redshifts (Az/(1 + z) < 2%). rest- 
frame colors and stellar population parameters. The 
SFRs presented in this paper are based in part on 
Spitzer-MIP S fluxes at 24/xm tha t are derived from the 
S-COSMOS (ISanders et al.lf2007h and FIDEL0 surveys. 
A comprehensive o vervie w of the survey can be found 
in I Wh itakcr e t al.l (|2011|) . The stel lar masses used i n 
this work are derived u sing FAST dKriek et al.l l27)M) . 
with iBruzual fc Chariot! (j2003!) models that assume a 
iChabrierl (|2003D initial mass function (IMF), solar metal- 
licity, exponentially declining star formation histories 
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and dust extinction following the iCalzetti et al.1 (j2000l ) 
extinction law. 

The SFRs are determined by adding the UV and IR 
emission, SFR X jy +m = 0.98 x 10- 10 (L m + 3.3L 2 soo) 
dKennicuttl I1998D adapted for the Kroupa IMF by 
iFranx et al.l (|200*8I ). accounting for the unobscured and 
obscured star formation, respectively We adopt a 
luminosity-independent conversion from the observed 
24/im flux to the total IR luminosity = L(8- 

1000/ini) ), based on a s i ngle t emplate that is the log av- 
erage of iDale fe Heloul (120021) tem plates with 1 < a < 
2.5, following IWuvts et al.l (120081): IFranx et al.1 (poM ): 
IMuzzin et al.l (|2O10h . and in good median ag reement with 
recent Herschel /PACS measurements by IWuvts et al.l 
(|2011al) . The luminosities at 2800 A (L2800) are derived 
directly from the best-fit template to the observed pho- 
tometry, u sing the same method ology as the rest-frame 
colors (see lBrammer et"al~ll2011l ). 

With accurate rest-frame colors, it is possible to 
isolate "clean" samples of star-forming and quies- 
cent galaxies using two rest-frame colors out to 
high redshifts tiLabbe et al.l [20051: IWuvts et al.l [20071: 
Williams etHI l2009Hllbert et al.1 120091 : iBrammer et al.1 
2011tlWhitake7 at al.ll2*011l ). The quiescent galaxies have 
strong Balmer/4000A breaks, characterized by red U-V 
colors and bluer V—J colors relative to dusty star-forming 
g alaxies at the same U- V color. 

iWhitaker et al.l (| 2 1 If) demonstrated that there is a 
clear delineation between star-forming and quiescent 
galaxies with the NMBS data set. Quiescent galaxies are 
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identified using the criteria U — V > 0.8 x (V — J) + 0.7, 
U—V > 1.3 and V — J < 1.5, and they are excluded from 
the bulk of this analysis. The sample of star-forming 
galaxies is selected independent of the SFR indicator and 
stellar population synthesis model parameters, enabling 
an unbiased measurement of the star formation sequence. 

3. THE STAR FORMATION SEQUENCE 

Figure Q] shows the star formation sequence, log(M*)- 
log(SFR), in five redshift bins out to z = 2.5. The 
greyscale represents the density of points for star-forming 
galaxies selected in Section [21 with the running median 
and biweight scatter color-coded by redshift. The mass 
completeness limits are estimated from the 90% point- 
source complet eness limits der i ved fr om the unmasked 
simulations by iWhitaker et al.l (|2011[) . The SFR com- 
pleteness limits correspond to the 3<r 24^tm detection 
limit (17.6 /xJy) at the highest redshift of each bin. All 
24^m detections < la are replaced with the la upper 
limit, resulting in a flattened tail of the log(SFR)-log(M*) 
relation at low M*, where the samples are incomplete. 

3.1. Quantifying the Star Formation Sequence 

The running medians and dispersions are measured for 
all star-forming galaxies, and those above the mass and 
SFR completeness limits are indicated with filled symbols 
in Figure [T] and fit with a power law, 

log(SFR) = aO)(logM* - 10.5) + f3(z) (1) 

We find that the slope a gradually evolves with red- 
shift, while the normalization f3 has a stronger evolution: 



a(z)=0.70-0.13z 
^(z)=0.38 + 1.14z-0.19z 2 



(2) 
(3) 



Previous studies have found hints of a similar trend 

for t he slope to flatten toward z ~ 1 or high M^ 

(e.g INoeske et al.ll2007l: iKarim et alll20lH IWuvts et all 
120121) . Although the slope is consistent with a gradual 
evolution toward shallower values, we note that the high- 
est redshift bins are only for galaxies with stellar masses 
> 10 10 7 Mq, due to incompleteness at lower masses. We 
will address this issue in Section 1X21 

The average SFR of star-forming galaxies steeply de- 
clines from z = 2.5 to the present day, changing by 
roughly 0.2 dex Gyr -1 , confirming the results of pr evious 
studies (e.g. INoeske et alj|2007t lElbaz et al.l 120071 etc.). 
The normalization /3 decreases by a factor of 4 from z — 2 
to z = 0, while the global SFR den sity has decreased a 
facto r of 6 over the same interval dHopkins fc Beacoml 
120061) . We note that j3{z) is the evolution in the specific 
SFR (sSFR ee SFR/M ) at log(MJ = 1 0.5 (consistent 
withlElbaz et al.ll2007t IDaddi et al.ll2008t IPannella et al.l 



[200a iDamen et al. l l2009HMagdis et alJl2010t etc). 

Interestingly, the scatter in the star formation sequence 
is observed to be roughly constant at ~ 0.34 dex, both 
with redshift and stellar mass, consistent with the idea 
that it reflects variations in the gas accretio n history. The 
observ ed scatter agrees with results from INoeske et al.l 
(|2007l) at z < 1 , although we note that there are hints of a 
lower scatter for the most massive galaxies in our highest 
redshift bin. We find that these results are robust against 
changes in sample selection. If we include all galaxies 
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Fig. 2. — More massive star-forming galaxies have higher 
log(Ljpt/Luv) ratios and redder rest-frame U—V colors, suggest- 
ing increasing amounts of dust attenuation. The dotted line is the 
mass-completeness limit. 



(quiescent and star-forming) with > la detections, the 
scatter increases by 0.07-0.11 dex. 

3.2. The Effects of Dust Attenuation 

Previous studies have noted that d ust attenuation is a 
strong function of stellar mass (e.g.. IReddv et al.l [2006L 
[20101 IPannella et al.l [20091: IWuvts et al.l 1201 lbt) . where 
the most massive galaxies are more highly obscured. We 
find a similar trend with the NMBS data in Figure [2j 
higher mass galaxies have larger Ltr/Ltjv ratios. The 
SFRs derived from L2800 alone exhibit a flat trend with 
stellar mass (or even hints of less UV flux toward higher 
M*), whereas SFRs derived from Lir show a strong trend 
with M*. The star formation-mass relation only mate- 
rializes when dust attenuation is properly taken into ac- 
count. 

From Figure [21 we also see that the average rest-frame 
U—V color of star-forming galaxies becomes redder with 
increasing stellar mass (and increasing Lir/Ltjv)- The 
red colors may be attributed to dust, although older stel- 
lar populations have similarly red colors. It is notoriously 
difficult to differentiate between subtle age and dust ef- 
fects with a single color alone. However, because Lir 
is independent of the rest-frame color measurements, we 
can isolate galaxies with red colors predominantly due to 
dust. 

The strong dependence of Lir/Ltjv and the rest-frame 
U-V color on stellar mass raises the question of whether 
the star formation sequence can be "resolved" into dis- 
tinct populations of star-forming galaxies. We first con- 
sider the following two samples at 1.0 < z < 1.5 (where 
we have a large number of galaxies at the high-mass end 
while maintaining a modest mass-completeness level): 1) 
blue star-forming galaxies with little dust obscuration 
(U -V < 1.3, log(Li R /L uv ) < 1.5), and 2) dusty red 
star-forming galaxies (U — V > 1.3, log(LiR/Luv) > 1-5). 

When we select blue star-forming galaxies, similar to 
iPeng et all (|2010D . the slope of the star formation sc- 
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Fig. 3. — The log(SFR,)-log(M*) relation for blue star-forming galaxies with low dust has a best-fit slope close to unity (left), whereas 
dusty, red star-forming galaxies follow a shallower slope (middle). Although selected independent of their measured SFRs and M*, blue, 
dusty galaxies and red, low-dust galaxies lie above and below the star formation sequence, respectively (right). X-ray detections are 
indicated with a black open circles, and those galaxies without X-ray coverage are indicated with crosses. 



quence is close to unity and the observed scatter de- 
creases to 0.25 dex (left panel of Figure [3]). However, 
dusty red star-forming galaxies are an increasing frac- 
tion of the massive ga laxy population toward z = 2 (e.g., 
iWhitaker etHl [2010 ) . The middle panel of Figure E2 
shows the star formation sequence for galaxies with red 
colors that we attribute to dust due to the high Lir/Ltjv 
ratios. This population of galaxies has a similarly small 
observed scatter of 0.25 dex, but a shallower slope of 
~ 0.8. 

While it has been shown t hat higher mass ga laxies tend 
to have lower sSFRs (e.g.. iZheng et al.l 120071 ). previous 
studies have not always distinguished between actively 
star-forming and passive stellar populations. Here, we 
find that among actively star-forming galaxies, the star 
formation sequence for red dusty (high-mass) galaxies 
has a shallower slope, as compared to that for blue low- 
dust (low-mass) galaxies. This result implies that both 
massive quiescent and star-forming galaxies have lower 
sSFRs and hence older ages. 

3.3. What causes galaxies to be outliers on the star 
formation sequence? 

We next consider the sample of star-forming galaxies 
with "anomolous" combinations of rest-frame color and 
Lir/Ltjv ratios to determine where these galaxies lie in 
the SFR-M* plane. We stress that we are not select- 
ing these galaxies to be outliers in the log(SFR)-log(M*) 
plane, rather to be unusual in their observed color and 
dust properties. The galaxies with red rest-frame colors 
and low Lir/Ltjv ratios fall in the lower envelope of the 
star formation sequence (right panel of Figure [3]) . These 
low sSFR galaxies may be in the process of shutting down 
star formation. Similarly, wc identify galaxies that have 
blue colors but high Ltr/Ltjv ratios. We find that dusty, 
blue galaxies have high sSFRs, predominantly residing in 
the upper envelope of the star formation sequence. This 
suggests that they are in a starburst phase and/or may 
have active galactic nuclei (AGN). 

To search for signatures of AGN, we match the NMBS 
catalogs within a 1 arcsecond radius to the X-ray point 
sourc e catalogs from the AEGIS-X survey (lLaird et al 



catalogs are indicated with open circles in Figure [3] We 
find that the fraction of galaxies with X-ray detections 
is ~ 4% in all populations of star-forming galaxies, im- 
plying that the central black holes may become activ e 
at any stage of star formation (e.g., (Saiitini e t al. 1 121)121 ). 
Interestingly, X-ray detected galaxies tend to be outliers 
in the log(SFR)-log(M*) plane (see Figure [3]), probably 
because a fraction of the 24/im emission is not associated 
with star formation, as assumed. 

4. DISCUSSION 

We find that the star formation sequence is not linear, 
with SFR oc M^' 6 and a constant observed scatter of 
0.34 dex. If we only select blue ga laxies, however, w e 
do find a linear relation, similar to iPeng et all (|2010( ) . 
This selection removes red, dusty star-forming galaxies 
at the hi gh mass end, which ha ve a shallower slope. We 
note that iPannella et all (|2009f ) also find a slope close to 
unity at z ~ 2, however this may be a result of the "BzK" 
selection. Here, we are able to analyze the star formation 
sequence with a mass-complete sample of galaxies. 

By accounting for the ratio of Lir/Ltjv (a proxy for 
dust) and the rest-frame U-V colors of star-forming 
galaxies, we develop a simple picture that describes how 
galaxies populate the log(SFR)-log(M*) plane. We have 
identified four distinct populations: 1) quiescent galax- 
i for log(IvL) > 10 at 1.0 < z < 1.5), 2) actively 



2009 ) and the Chandra COSMOS survey (jElvis et al 
20091) . Those galaxies with counterparts in the X-ray 



les 

star-forming galaxies with "normal" colors and associ- 
ated Lir/Ltjv ratios (54%), 3) red star-forming galaxies 
with low Lir/Ltjv ratios and low sSFRs (11%), and 4) 
blue star-forming galaxies with high Lir/Ltjv ratios and 
sSFRs (7%). 

Among the galaxies that populate the "normal" star 
formation sequence, we see a continuous sequence of in- 
creasing levels of dust attenuation with increasing stel- 
lar mass and a decreasing slope in the log(SFR)-log(M*) 
relation, i mplying decreasing s SFRs with stellar mass 
(Figure 0]). iWuvts et al.l (|2011bD observed similar trends 
across the SFR-M* plane, where Lir/Ltjv increases along 
the sequence and at a given M* towards high SFRs. The 
dependence of sSFR on M* appears to introduce a slight 
curvature to the star formation sequence. A power-law 
fit to all star-forming populations results in an observed 
scatter that is 0.09 dex larger than that of the "normal" 
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Fig. 4. — A cartoon (left) and observed (right) view of how star-forming galaxies populate the log(SFR)-log(M*) plane, with 54% of all 
galaxies residing on the "normal" star formation sequence, showing increasing amounts of dust (higher Ljr/Ltjv ratios) and lower sSFRs 
(shallower slope) towards higher stellar mass. 7% of galaxies have blue colors and high Ljr/Ltjv ratios, falling in the upper envelope (grey). 
11% of galaxies have low Lir/Ltjv ratios and red colors, populating the lower envelope (also grey). 



sequence, with a significantly shallower slope of ~ 0.6. 
We note that some models predict slopes that are too 
steep compared to the ob servations (e.g., IBouche et al.l 
120101 : IDutton et al.l 120101) ; it is possible that discrep- 
ancies between the inferred SFRs may be alleviated if 
the stellar IMF is systematically weighted toward more 
high-mass star formation i n rapidly star-forming galaxies 
(jNaravanan fc Pavel I20T2I ). 

The observed scatter of the "normal" star formation 
sequence is 0.25 dex, which includes contributions from 
both random and systematic errors. We estimate the av- 
erage random scatter by perturbing the 24/^m photome- 
try and photometric redshifts within the la error bars for 
100 realizations, finding a small contribution of ~ 0.05 
dex. Additionally, about 0.08 dex scatter is introduced 
because the average SFR evolves within the redshift bin. 
Finally, uncertainties in the conve rsion from 24/Ltm flux t o 
Lir introduce ~ 0.15 dex scatter ([MarciUac et al.lfeOOSl ). 
although we note th at this value is quite uncertain (e.g., 
iWuvts et al.ll2011al) . In total, we estimate that random 
and systematic errors introduce 0.18 dex scatter to the 
star formation sequence, from which we estimate the in- 
trinsic scatter to be 0.17 dex . Semi-analytic models pre- 
dict an even smaller scatter (0.12±0.1 dex). However, the 
model scatter may be underestimated due to a simplified 
treatment of the halo mass accretion history. Cosmolog - 
ical hydrodynamical simulations by iDekel et al.l ([2009D 
find a scatter of up to 0.3 dex in the gas accretion rates 
on to galaxies in 10 12 M© haloes at z = 2.5. If this scatter 
translates linearly into scatter in the SFRs, this may rec- 
oncile the differences between the models and observed 
scatter. 

To bolster the cartoon view presented in Figure |4l we 
additionally consider the composite SEDs of normal star- 
forming galaxies in bins of stellar mass. Due to the in- 
creasing levels of dust attenuation, we see a clear evo- 
lution of the composite SEDs in Figure [5] On average, 
the most massive star-forming galaxies have characteris- 
tically dusty spectral shapes, with a 2175A dust feature 
evident, whereas lower stellar mass galaxies have decreas- 
ing amounts of dust obscuration. We note that we see 



similar trends at higher and lower redshifts, but are un- 
able to make robust statements due to incompleteness. 

We speculate that the small fraction of dusty, blue 
galaxies are mainly starbursts, as very few appear to 
be associated with X-ray sources. Remarkably, we see 
that the spectral shapes of these galaxies are all very 
similar, irrespective of stellar mass (upper right panel in 
Figure [5]). This suggests that the same physical process 
is dominating the stellar light, possibly a merger driven 
starburst. 

While 28% of galaxies with log(M*) > 10 have already 
quenched their star formation at 1 < z < 1.5, we find 
that 11% may be in the process of shutting down star for- 
mation. These galaxies have red colors and low Lir/Ltjv 
ratios, occupying the lower envelope of the star forma- 
tion sequence. We compare the composite SEDs of these 
red, low-dust star-forming galaxies to that of quiescent 
galaxies at the same stellar mass and redshift (solid lines 
in bottom right panel of Figure [5]) . These galaxies have 
similar rest-frame V—J colors to quiescent galaxies, but 
somewhat bluer rest-frame U-V colors, consistent with 
the idea that they are in the process of shutting down star 
formation and may soon migrate to the red sequence. 

By studying four distinct populations of galaxies se- 
lected from the NMBS, we have demonstrated that quan- 
tifying the observed properties of the star formation se- 
quence and how the sequence evolves with time requires 
a thorough understanding of the selection techniques and 
biases. A consequence of the strong dependence of dust 
attenuation on stellar mass is that measurements of the 
star formation sequence will depend critically on the sam- 
ple selection. The gradual evolution we measure of the 
slope of the star formation sequence toward shallower 
values at high-z is driven by the combination of the cur- 
vature of the "normal" star formation sequence and the 
evolution of the mass-completeness limits with redshift. 
An in-depth analysis of the physical properties of these 
galaxies and comparisons between the observations and 
models will help constrain the physical mechanism driv- 
ing this potential curvature and the outliers from the star 
formation sequence. 
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Fig. 5. — The rest-frame composite SEDs of galaxies on the "normal" star formation sequence (left) showing increasing levels of dust 
attenuation with stellar mass. The spectral shape of dusty, blue galaxies appears to be independent of stellar mass (upper right). Galaxies 
in the process of shutting off their star formation (bottom right) show larger amounts of rest-frame UV emission relative to quiescent 
galaxies at the same stellar mass and redshift (solid lines). 
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